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By 
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May  1997 

Chairman:  Dr.  Joseph  Glover 
Major  Department:  Mathematics 

Since  Robert  Brown,  an  English  botanist,  observed  in  1828  the  highly  irregular 
motion  of  pollen  grains  suspended  in  water,  Brownian  motion  has  been  an  object 
of  study  in  both  pure  and  applied  mathematics.  From  a  probabilistic  standpoint, 
it  is  viewed  as  the  most  basic  stochastic  process  with  continuous  time  parameter. 
Its  application  is  involved  in  mathematical  models  of  random  phenomena  arising  in 
the  natural  environment.  From  an  analytic  standpoint,  Brownian  motion  is  used  to 
develop  powerful  tools  in  functional  analysis,  due  to  its  close  connection  with  Wiener 
measure  and  Newtonian  potential  theory. 

As  a  special  example  of  a  diffusion  process,  Brownian  motion  is  studied  for 
its  behavior  in  a  variety  of  state  spaces  with  particular  boundary  conditions.  Our 
interest  in  this  study  is  the  path-by-path  construction  of  a  diffusion  process  from 
Brownian  motion  in  Rd  where  d  >  3. 


v 


The  construction  is  based  on  inversions  with  respect  to  spheres  in  Rd.  We  in- 
vestigate certain  geometric  properties.  A  systematic  method  is  developed  to  partition 
the  state  space  into  a  series  of  closed  subsets.  A  superharmonic  function  h  is  defined 
on  Rd  which  has  an  invariance  property  under  inversions.  Brownian  motion  is  then 
transformed,  in  both  time  and  space,  while  its  fundamental  properties  as  a  diffusion 
is  retained. 

The  resulting  diffusion  process  has  normal  reflecting  boundaries.  In  the  interior 
of  its  state  space,  it  is  an  /i-Brownian  motion. 


CHAPTER  1 
INTRODUCTION 

1.1  Overview 

A  diffusion  process  is  a  strong  Markov  process  with  continuous  sample  paths. 
In  general,  a  diffusion  process  can  be  constructed  by  analytic  methods  or  probabilistic 
methods. 

The  analytic  methods  center  on  the  solution  of  a  parabolic  partial  differential 
equation  with  certain  coefficients.  The  major  difficulty  with  this  approach  is  that  it  is 
generally  a  non-trivial  problem  to  prove  that  this  parabolic  equation  has  a  sufficiently 
large  family  of  solutions  with  required  properties  to  characterize  the  associated  diffu- 
sion process.  Besides,  since  the  knowledge  about  this  process  can  only  be  obtained  by 
studying  the  associated  parabolic  equation,  there  is  a  loss  of  probabilistic  intuition. 

Another  approach  is  to  construct  the  diffusion  process  path-by-path,  on  the 
basis  of  an  underlying  probability  space  (0,  T,  P).  More  probabilistic  insights  can  be 
gained  by  this  approach  since  the  process  is  studied  directly  in  its  state  space.  One 
way  to  achieve  this  is  by  Ito's  method  of  stochastic  integral  equations.  Another  way 
is  to  construct  functions  of  well  studied  Markov  processes,  such  as  Brownian  motion. 
This  is  the  method  used  in  this  thesis. 

Specifically,  starting  from  Brownian  motion,  we  construct  a  multi-dimensional 
diffusion  process  which  is  an  /i-Brownian  motion  with  normal  reflecting  boundary. 
The  state  space  of  this  process  is  the  complement  of  n  separated  open  unit  balls  in 
Rd,  d  >  3. 

The  topic  of  Brownian  motion,  with  reflecting  boundary  has  been  studied 
with  both  probabilistic  and  analytic  approaches.  Results  on  one-dimensional  reflected 
Brownian  motion  can  be  found  in  Hida  [1]  and  Dynkin  [2].  In  Chung  and  Williams  [3], 
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they  discussed  the  difficulty  in  the  study  of  reflecting  Brownian  motion  in  two  or  more 
dimensions.  Analytical  theory  associated  with  these  processes  and  references  to  earlier 
articles  can  be  found  in  Harrison  and  Reiman  [4]. 

Glover  and  Rao  [5]  studied  reflected  Brownian  motion  in  R2  for  a  special  class 
of  domains  using  inversions.  In  that  article,  they  developed  a  systematic  partition 
of  R2.  A  Markov  process  is  then  constructed  path-by-path  though  time-changed 
Brownian  motion.  Finally,  this  process  is  identified  as  Brownian  motion  with  normal 
reflection  on  the  boundary  of  its  state  space. 

Their  work  is  the  motivation  of  this  study.  One  significant  difference  between 
R2  and  Rd,  d  >  3  is  that  Brownian  motion  is  recurrent  in  the  former  but  tran- 
sient in  the  latter.  This  difference  implies  that  no  reflected  Brownian  motion  can 
be  constructed  in  Rd,  d  >  3  by  an  appropriate  choice  of  time  change  and  inversions. 
However,  under  certain  constraints,  a  special  superharmonic  function  h  can  be  de- 
fined on  a  complement  of  a  polar  set  in  Rd  such  that,  for  an  appropriate  time  change, 
we  can  construct,  path-by-path,  an  /i-Brownian  motion  with  normal  reflection  on  its 
state  space. 

This  construction  is  the  main  result  of  this  thesis. 

1.2    Critical  Components 

The  critical  components  of  this  path-by-path  construction  pose  several  prob- 
lems that  have  not  yet  been  sufficiently  addressed  in  existing  literature.  These  prob- 
lems have  to  be  solved  in  order  to  complete  our  construction.  We  are  going  to  high- 
light these  critical  problems  in  this  section  and  start  rigorous  definitions  and  proofs 
in  Chapter  2. 

For  each  open  ball  Bj  =  B(cj,  =  1,2,...  ,n,  the  inversion  with  respect  to 
the  sphere  dBj  plays  an  important  part  in  the  construction.  It  is  defined  as 

Ij(x)  =  Cj  +  ||2(s  -  Cj) 
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for  x  €  Rd  -  {cj}.  Geometrically,  it  maps  a  point  x  G  Bj  to  a  point  lying  on  the  ray 
starting  at  Cj  and  passing  through  x,  such  that 

||/j(s)-q||x||s-c,||-l 

Obviously  Ij(x)  Bj.  If  x  &  Bj,  then  Ij(x)  G  and  the  same  identity  holds.  This 
collection  of  inversions  generates  a  non-Abelian  group  Q  of  homeomorphisms. 

The  first  problem  is  to  understand  the  properties  of  the  composition  of  these 
inversions.  Although  inversion  with  respect  to  a  single  sphere  has  been  introduced 
and  studied  in  the  literature,  the  geometric  properties  of  inversions  with  respect  to 
n  spheres,  where  n  >  1,  have  not  been  given  detailed  treatment.  It  is  necessary  to 
fully  investigate  these  properties,  which  are  not  always  trivial  due  to  the  non-Abelian 
nature  of  the  group  Q  . 

The  second  problem  is  also  associated  with  the  non-Abelian  nature  of  Q .  In 
R2,  Brownian  motion  is  recurrent  and  an  appropriate  time  change  will  result  in  a 
Markov  process  that  is  invariant  under  inversions  in  Q .  However,  since  Brownian 
motion  is  transient  in  Rd  where  d  >  3,  a  time  change  alone  is  no  longer  sufficient 
to  guarantee  this  invariant  property.  An  /z-transform  is  therefore  necessary  to  be 
coupled  with  the  time  change.  This  requires  that  the  superharmonic  function  h  has 
to  satisfy  the  identity 

h{x)  =  \\r -I  W-*h{Ij{x)) 
for  every  1  <  j  <  n.  Although  it  is  easy  to  find  a  superhamonic  function  that  satisfies 
this  condition  in  the  case  of  n  =  1,  the  existence  of  such  an  h,  even  in  the  case  of 
n  =  2,  is  not  obvious.  For  any  intuitive  choices  of  the  superharmonic  function,  the 
requirement  of  one  inversion  is  not  compatible  with  that  of  another.  The  difficulty 
here  is  to  find  one  superharmonic  function  h  that  simultaneously  satisfies  the  above 
identity  for  all  1  <  j  <  n. 

The  next  problem  involves  whether  a  many-to-one  continuous  function  maps 
a  Markov  process  to  a  Markov  process.  As  is  pointed  out  in  Glover  and  Mitro  [6],  this 


is  not  always  true.  The  particular  continuous  function  we  employ  has  to  be  studied 
to  insure  that  what  we  construct  is  indeed  a  Markov  process. 

Finally,  since  this  diffusion  is  constructed  path-by-path  from  a  Brownian  mo- 
tion, its  distribution  can  only  be  studied  through  that  of  the  Brownian  motion.  We 
have  to  characterize  the  semigroup  of  this  diffusion  by  computing  its  infinitesimal 
generator  based  on  its  connection  with  the  original  Brownian  motion. 

1.3    Outline  of  the  Construction 

Chapter  2  covers  the  study  of  the  geometric  properties  of  inversions  and  the 
state  space.  Definitions  and  basic  properties  are  discussed  in  Section  2.1.  In  Section 
2.2,  properties  of  Q  are  investigated.  The  equivalence  relationship  introduced  by  Q  and 
a  systematic  partition  of  Rd  is  studied  in  Section  2.3. 

Chapter  3  is  devoted  to  the  construction  of  two  functions.  In  Section  3.1,  we 
define  a  positive  function  g  which  introduces  a  time  change  later.  Another  function  h 
is  defined  in  Section  3.2  and  its  superharmonicity  is  established  in  Section  3.3.  Other 
properties  of  this  function  are  also  demonstrated  in  these  two  sections. 

Chapter  4  contains  the  main  result  of  this  thesis.  We  give  the  path-by-path 
construction  of  a  diffusion  Xt  in  Section  4.1  and  then  identify  the  infinitesimal  gen- 
erator of  its  semigroup  in  Section  4.2.  The  invariance  and  symmetry  properties  of 
this  process  are  studied  in  Section  4.3  and  4.4.  Finally,  in  Section  4.5,  a  Markov 
process  Wt  is  constructed  path- by-path  through  Xt.  The  infinitesimal  generator  of 
the  semigroup  of  Wt  is  computed.  The  result  demonstrates  that  Wt  is  an  /i-Brownian 
motion  in  the  interior  of  its  state  space  and  is  normally  reflected  on  the  boundary. 

Several  restrictions  are  imposed  in  the  process  of  this  construction.  In  Section 
4.6,  the  possibility  and  difficulty  of  relaxing  them  are  discussed. 


CHAPTER  2 
THE  STRUCTURE  OF  THE  STATE  SPACE 

2.1    Definitions  and  Basic  Properties 

In  this  section,  we  shall  define  a  family  of  inversions  on  Rd  ,  the  one  point 
compactification  of  d-dimensional  real  space.  Each  of  these  inversions  is  a  homeo- 
morphism  of  the  space  onto  itself. 

For  any  open  ball  B  C  Kd,  let  dB  be  the  boundary  of  B.  From  now  on,  unless 
otherwise  specified,  we  shall  refer  to  dB  as  the  sphere  of  B. 

Let  us  denote  by  Rd  the  one-point  compactification  of  Rd: 

W  =  Rd  U  {00} 

For  any  point  c  G  Rd,  let  B  =  B(c,  1)  be  an  open  ball  of  radius  1,  centered  at  c.  We 
define  a  map  /:  Rd  ->  Rd  as  follows. 
Definition  2.1.1.  I :  Rd  ->  Rd  is  defined  as 

lit)   =  c+yr^ — "7|o  for:r£{c,oo} 
1(c)   =  00 
J(oo)   =  c 

This  map  leaves  each  point  on  the  sphere  of  B  fixed  while  it  maps  each  interior 
point  of  B  into  B°,  and  vice  versa;  i.e.,  if  x  G  iT,  then  I(x)  G  B.  Thus,  /  is  an 
inversion  with  respect  to  the  sphere  of  B. 

From  this  definition,  we  can  derive  the  following  basic  properties  of  /. 


Proposition  2.1.1.      1.  /  is  idempotent,  i.e., 

r1  =  /  (2.1) 

2.  For  all  x  g  {c,  oo}, 

\\x-c\\x\\I(x)-c\\  =  l  (2.2) 
Proof.  These  are  direct  consequences  of  the  definition. 
1.  By  definition,  for  all  x    {c,  oo} 

/(/(*))  =  /(e  +  ^fi) 


x-c  Hi-cll2 
C+(llx-clPMl|x-cl|*' 


=  X 


Moreover,  we  have 

/(/(c))  =  /(oo)  =  c 

and 

/(/(oo))  =  /(c)  =  oo 

So,  (2.1)  is  true. 

2.  For  any  x  0  {c,  oo},  0  <  ||x  -  c\\  <  oo,  0  <  \\I(x)  -  c||  <  oo,  and 

||x  -  c\\  x         -  c\\   =   \\x  -  c\\  x  \\c  +    *_  ^    -  c\ 

ii        ii  \\x-c\\ 


x  —  ell2 


=  1 

So,  (2.2)  is  true. 


□ 


In  order  to  study  diffusion  processes  later  using  the  inversion  /,  we  need  to 
learn  more  about  its  properties.  To  do  this,  we  must  first  prove  the  following  two 
lemmas  in  geometry.  Properties  of  /  are  studied  in  Propositions  2.1.2  and  2.1.3  and 
in  the  next  section. 

Lemma  2.1.1.  Suppose  two  straight  lines  passing  through  points  A,  D  and  B,  C 
respectively,  intersect  at  point  O.  See  Figure  1. 

D 

A  \ 


u  B 

Figure  1 

If 

OBOC  =  OA-  OD 

then 

AB  =  OA  =  OB 

CD     OC     OD  [  '  ' 

Proof.  Consider  AOAB  and  AOCD.  By  the  given  condition,  we  have 

OA  OB 


OC  OD 

Since  ZAOB  =  /.COD,  we  can  conclude  that  AOAB  and  AOCD  are  similar.  Hence 

AB     OA  OB 


CD     OC  OD 
Obviously,  this  is  also  true  in  the  special  case  when  OA  =  OD. 


Lemma  2.1.2.  Suppose  two  straight  lines  passing  though  points  A,  D  and  B,  C 
respectively,  intersect  at  point  O  (See  Figure  2). 


If 


then 


Figure  2 


OB  OC  =  OA-  OD 


AC     OA  OC 


(2.4) 


BD  ~  OB  ~~  OD 

Proof.  Consider  AOAC  and  AOBD.  By  the  given  condition,  we  have 

OA  _  OC 
OB  ~  OD 

Since  ZAOC  =  ZBOD,  we  can  conclude  that  AOAC  and  AOBD  are  similar.  Hence 

AC  _OA  _  OC 
BD  ~  OB  ~  OD 


Obviously,  this  is  also  true  in  the  special  case  when  OA  =  OD. 


□ 


Equipped  with  Lemmas  2.1.1  and  2.1.2,  we  are  ready  to  study  the  following 
properties  of  the  inversion  /. 
Proposition  2.1.2.  For  all  x,  y  0  {c,  oo}, 

Ik  -  y\\ 


\\I(x)-I(y)\\ 


\x-c\\  x  \\y-c\ 


(2.5) 
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Proof.  Depending  on  the  relative  positions  of  x  and  y,  the  proof  can  be  divided 
into  the  following  cases. 

1.  If  both  x  and  y  are  on  the  sphere  of  B,  we  have  I(x)  —  x,  I(y)  —  y  and 
\\x  —  c||  =  \\y  —  c||  =  1.  The  proof  is  trivial. 

2.  If  both  x,  y  €  B  and  at  least  one  of  them  is  not  on  the  sphere,  we  have  two 
triangles  formed  by  points  {c,x,y}  and  {c,I(x),I(y)},  respectively,  as  shown 
in  Figure  3. 


Figure  3 

By  (2.2)  \\I(x)  -  c\\  x  ||x  -  c\\  =  1  =  \\I{y)  -  c||  x  \\y  -  c\\.  It  follows  from 
Lemma  2.1.1  that  these  two  triangles  are  similar;  hence 

II*  -2/11  \\y-c\\ 


or,  equivalently, 


Since 


we  get 


Therefore, 


||/(x)-/(y)||  ||/(x)-c|| 

P(*)-/(y)ll  _  \\m-4 


\x-y\ 


-  c||  = 


x  -  c 


\\I(x)  -  I(y)\\ 


\x-y\ 


\x-c\\  x  y-c 


wm-mw 


f  -  y\ 


\x  -  c||  x  \\y  -  c| 
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3.  If  both  x,y  £  B°,  and  at  least  on  of  them  is  not  on  the  sphere,  then  the  two 
triangles  formed  by  {c,x,y}  and  {c,  I(x),  I(y)},  respectively,  are  as  shown  in 
Figure  4. 


Figure  4 


In  this  case,  we  also  have  ||  J(a;)  -  c||  x  |}x  -  c\\  =  1  =  \\I(y)  -  c\\  x  \\y  -  c\ 
Lemma  2.1.1  can  be  applied  here  and  we  have 

||/(*)-/(y)||  _  \\m  -  4 


\x-y\ 


By  similar  argument,  we  get 


\\I(x)-I(y)\\  = 


\y-c\ 


\x  -  y\ 


\x-c\\  x  \\y-c\ 


4.  If  only  one  of  the  points  is  in  B,  then  the  triangles  formed  by  {c,  x,  y}  and 
{c,  I(x),I(y)},  respectively,  are  as  shown  in  Figure  5. 


Figure  5 


11 


In  either  case,  \\I{x)  -  c\\  x  \\x  -  c\\  =  1  =  ||/(y)  -  cj|  x  ||y  -  c||,  by  (2.2). 
According  to  Lemma  2.1.2,  we  have 

l|/(*)-/(y)ll  =  \\Hx)-4 
\\x-y\\  \\y-c\\ 

Hence,  by  the  same  argument  as  above,  we  have 

\\x-y\\ 


\\I(x)-I(y)\\  = 


\x-c\\  x  ta-c 


□ 


Proposition  2.1.3.  For  any  x,  y  <£  {c,  oo}, 

\\i(x)-y\\  =  \\y-4  f2  6) 

II* -ell 

Proof.  There  are  four  possible  cases  defined  by  the  relative  positions  of  x  and  y. 
They  are: 

(a)  x  e  B  and  y  &  B; 

(b)  x  &  B  and  y  e  B; 

(c)  x  €  B  and  y  €  £; 

(d)  x  g  B  and  y  g  B; 

as  shown  in  Figure  6. 

In  each  case,  we  have 

-  c||  x  ||y  -  c||  =  1  =  \\I{x)  -  c\\  x  \\x  -  c\\ 

Applying  Lemma  2.1.1  in  case  (a)  and  (b),  or  Lemma  2.1.2  in  cases  (c)  and  (d),  we 
have 

-y\\  _  \\y-4 


\x-I(y)\\  \\x-c\ 


(b) 

Figure  6 

Therefore,  we  have  the  identity  (2.6). 


(d) 


□ 


Proposition  2.1.4.  /  is  a  continuous  map  from  Rd  to  Rd  . 
Proof.  Let  us  verify  directly  that  /  is  continuous  at  c  and  oo. 

First,  suppose  a  sequence  {xn}™=l  in  Rd  has  point  oo  as  its  limit  point.  Then 
||x„||  -»  oo  as  n  ->  oo.  By  (2.2),  for  any  n, 

\\I(xn)  -  /(oo)||  =  ||/(a:„)-  cfl 

1 


< 


l 


\Xr 


Since  ||arn||  ->  oo  as  n  -¥  oo,  the  right  hand  side  approaches  zero.  Hence,  /  is 
continuous  at  oo. 
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Next,  suppose  a  sequence  {xn}™=l  in  Rd  converges  to  point  c,  i.e.,  ||rrn-c||  -*  0 
asn  -*  oo.  Then,  ||/(xn)||  +  ||c||  >  ||/(x„)-c||  =  H^-cJI-1  oo.  Since  ||c||  is  finite, 
we  have  ||/(a;n)||  ->  oo,  i.e.,  I(xn)  ->  7(c).  Therefore,  I  is  also  continuous  at  point  c. 

Generally,  let  {a^}^  be  a  sequence  in  Rd  that  converges  to  x0  ^  c.  Then 
||x0-c||  >  0  and  there  exists  e  >  0  such  that  ||xn-c||  >  e  for  all  n.  Since  ||xn-a;o||  ~*  0 
as  n  -»  oo,  by  (2.2), 

iir/    \     r/    mi  llgn-go||         ^  IK-Zoll  n 

"  /(Xo)l1  =  IK-cllxIbo-cH  -  ° 


as  n  ->•  oo.  So,  I  is  continuous. 


□ 


The  definition  of  I  and  Proposition  2.1.1  guarantee  that  I  is  a  bijection.  Thus, 
7  is  a  homeomorphism  of  Rd  onto  itself. 

Generally,  for  each  Bj  =  B(Cj,  1),  j  =  1,  2,  n,  we  can  define  an  inversion  Ij 
with  respect  to  dBj,  and  all  of  the  properties  discussed  above  hold  for  each  inversion, 
respectively. 

2.2    The  Group  Q  Generated  By  Inversions 

In  this  section,  a  finite  number  of  inversions  are  defined  in  the  same  way  as  in 
section  2.1,  and  the  group  of  homeomorphisms  of  Rd  generated  by  these  inversions  is 
discussed. 

Let  Bj  =  B(cj,  1),  3 ia  1,  2,  . . .  ,  n,  be  separated  open  balls  in  Rd.  By  this, 
we  mean  B~{  n  B~j  =  0,  for  t  ^  j.  Let  Ij  be  the  inversion  with  respect  to  dBj,  the 
sphere  of  Bj,  defined  in  the  same  way  as  7  in  section  2.1,  for  j  =  1,  2,  . . .  ,  n. 

Each  inversion  Ij  is  a  homeomorphism  of  Rd  onto  itself,  so  they  generate  a 
non-Abelian  group  Q  of  homeomorphisms  on  Rd  . 

Next,  we  shall  prove  a  fact  concerning  the  identity  of  this  group. 
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Lemma  2.2.1.  Let 

where  k  >  1  and  rij  G  {1,  2,  . . .  ,  n}  for  1  <  j  <  k.  If  n,  /  ni+1  for  1  <  j  <  k  -  1, 
then  S/e. 

Proof.  By  the  definition  of  /j,  for  any  x  &  Bj,we  have  Jj(x)  G  Bj. 

Suppose  5  =  Jni  J„2 . . .  Ink_Jnk  and     ^  ni+1  for  1  <  j  <  k  -  1.  Let  z  =  00, 
xfc  =  Jfc(z)  and 

Xj  =  ityfo+i) 

for  j  =  k  -  1, k  -  2,. . .  ,3,2,1.  It  follows  that  xk  =  Inis(oo)  G  Since  nk  ^  nfc_i, 
xfc  £  5^17.  Therefore,  zfc_i  =  Inh.iM  G  Repeating  the  same  argument,  we 

have 

for  all  1  <  j  <  k  -  1.  Thus,  xx  =  S(z)  G  Bm,  while  00  g"  JBni,  and  we  conclude  that 

□ 

Since  Q  is  generated  by  {A,  I2, . . .  ,  4},  any  5  G  £  -  {e},  can  be  represented 
as  a  finite  product  of  these  inversions. 

By  the  idempotent  property  of  the  inversions,  for  any  S  G  Q ,  if  one  represen- 
tation of  S  is 

we  can  always  find  an  equivalent  representation 

S  =  .  .  .  In' k, _1In' k, 

such  that  n'i  ^  n'i+i  for  every  1  <  i  <  k1  —  1. 

Next  we  shall  prove  that  this  representation  is  unique. 
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Lemma  2.2.2.  If  S  6  Q  -  {e}  has  a  representation  S  —  InJn2  ■  ■  ■  hk-Jnk  for  some 
ni,  n>2,  •  •  •  >  "A  and  nt  ni+i  f°r  every  1  <  i  <  k  —  1,  then  this  representation  is 
unique. 

PROOF.  This  lemma  can  be  proved  by  way  of  contradiction. 

Suppose  this  is  not  true,  and  we  have  another  representation 

where  m{  ^  mi+i  for  every  \  <  i  <  p  -  1  and  n{  /  ni+i  for  every  1  <  i  <  k  -  1. 
Without  loss  of  generality,  we  may  also  assume  that  p  <  k. 

We  claim  Ini  =  Jmi,  or  else  we  shall  have  a  contradiction. 

Since  5  =  ImiIm2  ■  ■  ■  Imp-Jmp,  and  every  inversion  is  idempotent,  it  follows 

that 

S      =  Imp  ^rrip- 1  ■  ■  •  Im2Im\ 

If  Jni  ^  Jmi ,  then  the  following  product 

T  =  ImpImp-\  ■  ■  ■  An2 -^mi^ni-^  •  •  ■  Ink-ilnk 

has  the  property  that  no  two  consecutive  inversions  are  the  same.  Therefore,  by 
Lemma  2.2.1,  T^e.  On  the  other  hand,  T  =  S^S  =  e.  This  contradiction  shows 
that  Ini  =  I  mi  ■ 

Now  we  have 

Imi  ■  ■  •  Imp-\Imp  ~  A12  ■  •  ■  Ink-\Ink 


By  mathematical  induction, 


Ifii  —  IfUi 


for  all  1  <  i  <  p. 

If  p  =  k,  the  proof  is  complete.  If  p  <  k,  we  have 

e  =  Inp+i  ■  •  ■  Ink 

But  since  ^  ni+i  for  every  1  <  t  <  k  —  1,  this  is  impossible  according  to 
Lemma  2.2.1. 

Hence,  the  representation  of  any  S  e  Q  -  {e}  is  unique. 

□ 
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From  now  on,  unless  otherwise  specified,  we  always  refer  to  this  unique  repre- 
sentation whenever  S  is  mentioned.  Now  that  we  have  the  unique  representation  of 
any  element  of  Q  —  {e}  ,  we  can  define  its  length. 

Definition  2.2.1.  For  any  S  =  InJn2 . .  .Ink_Jnk  €  Q  -  {e},  the  length  L  of  S 
is  defined  as  the  number  of  inversions  in  its  unique  representation.  L(S)  =  k  for 
S  =  JBlJfta . . .  Ink-Jnk-  K  we  define  L(e)  =  0,  then  L  :  Q  -¥%. 

The  following  is  a  generalization  of  Proposition  2.1.3. 
Proposition  2.2.1.  Let  S  e  Q  -  {e}.  For  any  x,  y  &  {ci,c2,...  ,cn,oo}, 

\\S(x)-y\\  _  \\y  -  S{oo)\\  (2  J) 


\\x-S-'(y)\\  \\x-S-^)\\ 
PROOF.  Suppose  S  =  IniIn2 . . .  Ink-Jnk  is  its  unique  representation.  The  conclusion 
can  be  proved  by  mathematical  induction  on  the  length  of  S. 

(a)  In  the  case  of  S  =  Jni,  it  follows  that  S'1  =  Ini  and  S(oo)  =  5_1(oo)  =  c^, 
and  Proposition  2.1.3  can  be  applied  directly  to  get  the  desired  identity. 

(b)  Assume  the  identity  is  true  for  Si  of  length  k  -  1;  i.e.,  for  Si  = 
we  have 

\\SM-v\\  fe-frMH 


\\x-s;\y)\\  IN-^MH 

Let  S  =  SJn.,  where  nk  ±  nk-V  Then,  S'1  =  InkSi  \  and,  by  the  inductive 


assumption, 


\\S{x)-y\\  =  \\SJnk(x)  -  y\\ 

=  l»  -  s'(00)l1  x  !<„(.)  -  sr 'Mil 


It  follows  from  (2.6)  that 

\\Ink(x)-S^(y)\\  =  \\x-InkS^(y)\\x 


|x  -  /nfc(00 
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and 


i/„l(l)  -  sr'(°°)ll  =  II*  -  /^sr'Mii  x  "^ffi 


Hence,  the  ratio 


\\Ink(x)  -  S{\y)\\  \\x-S-Hy)\\x\\S^(y)-Ink(oo)\\ 
\\Ink(x)-S;l(oo)\\     ||x-5-1(oo)||  x  \\S;\oo)-Ink(cx>) 


It  then  follows  that 

HO/x    „     „    C/  miv  II*  -^(y)!!  x  Mggll! 

H5(x)-y||   =  5,(00)11  x||x_5-1(oo)||X||5ri(oo)_^(oo)„ 

Ik  ~  S-l(y)\\      K\y)  -  Ink(oo)\\  x  \\y  -  5,(00)11 
||ar-  5-l(oo)||  \\S?(oo)  -  /Bfc(oo)|| 

Applying  the  assumption  again  for  5f\  we  have 

ll^-^MII  =  11^-^(00)11 

I,      c^n  v  11^(00) -^(oo)ll 
■  l|y_  5(00)11  X      ||y- 5,(00)11 

So, 

l^rHy)  -  /nt(oo)||  x  ||y  -  5,(00)11  -  ||y  -  5(oo)||  x  ^(00)  -  J„4(oo) 
Thus, 

nct\      ll         l|x-5-1(y)||      ||y  -  5(oo)H  x  \\S^(oo)  ~  ^Mll 
2/11   ~    Har- 5-i(oo)||  H^oo)  -  IBJk(oo)|| 

-  H*    5    (y)||x  |)a;_5_1(oo)|| 
Therefore,  we  have  (2.7). 


So,  by  mathematical  induction,  identity  (2.7)  is  true  for  5  of  any  finite  length. 

□ 


18 


Similarly,  Proposition  2.1.2  can  also  be  generalized  for  general  S  €  Q  —  {e}. 
In  doing  so,  for  any  S  G  Q  —  {e},  we  introduce  a  factor  Ds-  This  factor  is  useful 
in  our  understanding  of  the  properties  of  the  group  Q ;  it  is  also  important  in  the 
construction  of  a  superharmonic  function  h  that  will  be  discussed  in  the  following 
chapter. 

Proposition  2.2.2.  Let  S  =  InJn2 .  ..Ink_Jnk.  For  any  x,y  g  {c^,  c„2,  . . .  ,  c^,  co}, 
the  following  identity  is  true. 

™  ~ SWI1  =  ii,-5-(cc)ii!to-5-(0o)iixPg»  (2-8) 

where,  for  k  =  1, 

Ds  =  l;  (2.9) 

and  for  k  >  2, 

Ds   =   ||/ni(oo)-/„2(oo).../n,(oo)||  x  ||/n2(oo) -/n3(oo).../„fc(oo)||  x  •••  x 

||/„fc_a(oo)  -  V.M/^MH  x  HV^oo)  -  /B4(oo)||  (2.10) 

Proof.  This  result  can  be  proved  by  mathematical  induction  on  the  length  of  S. 

(a)  In  the  case  of  k  =  1,  S  =  Ij,  j  €  {1,  2,      n},  we  have  Proposition  2.1.2,  i.e., 

(b)  In  the  case  of  k  =  2,  suppose  5  =  IJj,  where  i,  j  G  {1,  2,      n}.  It  follows  that 
S-1  =       and  Ds  =  ||/<(oo)  -  /j(oo)||. 

By  Proposition  2.1.2,  we  have 

||5(x)-5(y)||   =   Wl^x)  -  UjMW 


\\Ijix)-  /i(oo)||  x  HJjfo)  -  J,(oo) 


Using  Proposition  2.1.2  again,  we  observe 
Moreover,  by  Proposition  2.1.3 

Pi(oo)  -7,(oo) 


\\Ijix)  -/i(oo)||    =    ||x- 1^(00)11  X      ||x_/  (oo)|| 

Combining  these  two  expressions,  we  have 

||/,(x)-/,(oo)||x  \\I,(y)  -/i(oo)|| 

=  |  -  WMi  x  J,  -  y.Mfl  x  n.Jg^-y^; 

Therefore, 

l|5(x)-5(y)"   =    Hx-^JJx  J- 7,(00)11  X 

Us  -  4(00)11*110-7,(00)11 


\\x  -  7^(00)11  x  ||y  -  7^(00)11  x  117,(00)  -  7,(oo)||2 

 \\x-y\\  

||ar  -  7,7,(00)11  x  ||y  -  7,7,(00)11  x  p^oo)  -  7J-(oo)||* 


||ar- S-^oo)!!  x  \\y-S-i(oo)\\xDs2 
(c)  Assuming  that  identity  (2.8)  is  true  for  S  of  length  fc— 1,  i.e.,  for  Sr  —  7n2 . . .  In 
we  have 

\\x  -  y\\ 

l|SlW  "  s,m  =  ik-sr'Miixii.-sr'Miixuj, 

Let  S  =  7„,Si  =  7ni7„2 . . .  Inic,  where  nx±n2.  Then  we  have 

S    =  iSj  7ni 

and,  by  (2.10) 

7>s   =   ||7ni(oo)-7n2(oo)...7nt(oo)||  x  ||7na(oo)  -  7„s(oo) . .  .7nt(oo)||  x  ••• 
||7„t_2(oo)  -  7^^(00)7^(00) ||  x  117^(00)  -  7„fc(oo)|| 
=  ||7Bl-5i(oo)||xDSl 


Now,  by  Proposition  2.1.2 


\\S(x)  -  S(y)\\   =   \\IniSi{x)  -  InMv)\\ 

\\Si(x)  -  S1(y)\\ 


||5i(x)-/ni(oo)||x||51(y)-/Bl(c») 


and  by  the  inductive  assumption, 

\\x  -  y\\ 

"Sl(l)  "  S'iM  =  ||x-Sr'(oo)!|x||K-Sr'(oo)||xD^ 
By  Proposition  2.2.1 

I4M-A.MI  =  iN-gr1/..(°°)llxliys:1gy 

l|Si(»)-/«M||  =  ||»-Sr1/.,(°o)||x'ly^~) 

Substituting  the  above  expressions,  we  obtain 
l|5(x)-5(y)||  = 


II*  -y\\ 

\\x- 

\oo)\\  x\\y-S^ (oo)\ 

1  x  DSl2 

llx-^MH 

\\x  — 

$i 

^(ooJIjxpn^OO)- 

Si(oo)\\ 

Ijy-^MH 

\\y- 

St 

1/»1(oo)||x||41(oo)- 

S1(oo)\\ 

\\x  - 

-y\\ 

\\x- 

s- 

■x(«))||  x  ||y-5r-i(c») 

II  x  DSl2 

Ik -2/11 

\\x- 

s- 

i(cx))||x||y-5-i(oo) 

II  x  DS> 

X 


Therefore,  identity  (2.8)  is  true  for  all  S  €  Q . 
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2.3    Orbits  and  the  Horizon  Set 

As  we  have  shown  in  section  2.1,  any  single  map  Ij  defines  an  inversion  in  Rd  . 
We  are  now  ready  to  investigate  the  geometry  of  these  inversions,  particularly,  the 
geometric  characteristics  of  elements  of  the  group  Q .  The  following  construction  is 
based  on  the  same  method  first  given  in  [5]  by  Glover  and  Rao  in  the  2-dimensional 
case. 

The  group  Q  of  homeomorphisms  on  Rd  defines  an  equivalence  relationship  in 
Rd  in  the  following  way: 

x  ~  y  if  and  only  if  y  —  S(x)  for  some  S  G  Q 

An  equivalence  class  is  formally  defined  as  an  orbit  below. 

Definition  2.3.1.  For  any  x  G  Rd  ,  the  orbit  of  x,  denoted  by  [x],  is  defined  as  the 
following: 

[x]  =  {S(x):  SeQ) 

Since  each  element  of  the  group  Q  is  a  homeomorphism  on  Rd  ,  every  point  in 
Rd  belongs  to  a  certain  orbit.  Proposition  2.3.1  below  characterizes  the  collection  of 
all  orbits.  First,  we  need  a  lemma. 

Lemma  2.3.1.  For  every  x  G  R1  -  D]=lBj,  and  for  every  S  eQ  such  that  L{S)  >  2, 
we  have 

S(x)  €  ^=1Bj 

Proof.  For  any  j  e  {1, 2, . . .  ,n},  let  Ij  be  the  inversion  defined  in  Definition  2.1.1 
with  respect  to  OBj,  where  Bj  =  B(cj,  1). 

By  definition,  if  x  G  dBj,  then  Ij(x)  =  x;  otherwise,  for  any  x  &  Bj, 
I3(x)  e  Bj. 

First,  suppose  S  =  Ink-Jnk-  For  every  x  e  Rd  -  U"=1Bj,  there  are  two 
possibilities.  If  x  ^  dBnk,  then  Ink{x)  G  Bnk.  Otherwise,  if  x  G  dBnk,  then 
Ink(x)  =  x  G  dBnk. 


But  since  rik-i  /  n*,  by  our  assumption, 

Bnk_,  n  B^k  =  0 

Therefore,  in  either  case,  Ink{x)  &  Bnk_v  So,  Ink_Jnk(x)  G  -B„fc  l. 

Now  we  can  prove  the  general  case  by  mathematical  induction.  Assume  for 
any  Sk  =  InJn2  ■  ■  -Ink.Jnk  with  2  <  k  <  N,  it  is  true  that  Sk(x)  G  Bni,  for  any 
iftf  -  Unj=lBj.  Then  for  5  =  InJnJn2  ■  ■  ■  InN-JnN,  where  n0  ^  m,  let 

we  have,  by  our  assumption,  y  G  Bni.  Hence,  y  £  £?„„.  It  follows  that 

S{x)=Ino(y)eBnoCU]=lBj 

By  mathematical  induction,  for  every  S  G  Q ,  and  every  xetf  -  U?=1#j,  we  have 
,S(:r)  G  UJUB,. 

□ 

This  lemma  leads  to  the  following  conclusion. 
Proposition  2.3.1.  Each  orbit  contains  at  most  one  point  in  Rd  -  Uj=lBj. 
Proof.  Suppose  both  x  and  y  are  in  Rd  -  U"=1Bj,  and  they  are  both  in  the  same 
orbit,  i.e.,  [x]  =  [y].  Then  there  exists  S  G  Q  with  S(x)  =  y,  We  can  show  that 
x  =  y. 

We  discuss  all  possible  cases  according  to  the  length  of  S. 
First,  by  Lemma  2.3.1,  the  case  of  L(S)  >  2  is  impossible  since 

xeW  -u]=lBj 

would  imply 

y  =  S(x)  G  U^Bj 


23 


Second,  let  us  assume  that  L(S)  —  1,  so  S  =  h,  for  some  k  G  {1,  2,  . . .  ,  n). 
If  x    dBk,  then  y  =  Ik{x)  G  .Bfc,  but  this  contradicts  the  assumption  that 

yeW  -  U]=1Bj 

So  we  must  have  x  G  dBk.  It  then  follows  that 

y  -  S{x)  =  Ik(x)  =  a; 

Finally,  if  L{S)  =  0,  then  S  —  e,  so  x  —  y.  □ 

In  general,  there  are  two  types  of  orbits.  The  first  type  contains  a  point  in 
Rd  -  U"=1Bj  and  the  second  type  does  not  contain  such  a  point.  Letting 

E  =  W  -U"=iBj 

we  can  define  the  horizon  set. 

Definition  2.3.2.  The  horizon  H  is  defined  as 

H  =  {x  <E  W  :[x]r\E  =  ®} 

Restricting  each  homeomorphism  of  Q  to  Hc,  we  can  see  that  the  group  Q  defines 
an  equivalence  relation  on  Hc  and  that  each  equivalence  class  contains  exactly  one 
point  in  E.  Let  Hc/Q  be  the  collection  of  all  equivalence  classes,  or  orbits,  in  Hc. 
For  any  [x]  in  Hc/Q ,  let  y  be  this  unique  point  in  E.  If  we  define  </> :  H°/Q  ->  E  by 

4>([x])  =  y 

then  <f>  is  a  bijection  from  H°/Q  onto  E.  Moreover,  if  we  define  a  map  $  :  Hc  — ►  E 

by 

=  <f>([x])  =  y 

then  $  is  a  map  from  Hc  onto  E. 
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Next  we  proceed  to  partition  Rd  into  the  union  of  a  sequence  of  regions.  Define 

Es   =   {S(x)  :  x  G  E}  for  any  S  G  Q 
Gk   =   {S  6  G  :  L(S)  =  k} 
Fk  =  UsegkEs 

Since  E  is  a  closed  set  and  every  S  G  G  is  continuous,  2?s  is  also  closed. 

For  k  —  0,  Gk  —  {e};  and  for  every  k  >  1,  the  cardinality  of  &  is 
n(n  -  l)*-1.  Therefore,      is  a  closed  set  since  it  is  a  finite  union  of  closed  sets. 

The  following  lemma  further  characterizes  the  horizon  set  and  its  complement. 

Lemma  2.3.2.  H  —  nSeg  S(\Jj=1Bj)  is  compact  and 

H°  =  Useg  ES 

Proof.  It  is  better  to  prove  the  result  for  Hc  first.  By  the  definition  of  H, 

Hc  =  {xeW  :[x]n£/D} 

Therefore,  for  any  S  G  G ,  and  any  x  G  Es,  we  have  S1-1^)  G  E  by  the  definition 
of  Es.  So,  [  x  ]  HE  ^  0.  Thus,  x  G  Hc  for  all  x  G  Es.  Therefore,  Es  C  Hc  for  all 
S  GG  ■  This  implies  that 

US€gEsCHc 

On  the  other  hand,  for  any  x  G  Hc,  [x]  n  E  ^  0.  So,  there  exists  5  G  £  such 
that  y  =  S(x)  G  E.  Since  a;  =  5_1(y)  and  y  e  E,we  know  that  a;  G  £"5-1.  Therefore, 

#c  c  Useg  Es 

Thus  we  have  the  desired  result. 

Next,  let  us  examine  the  set  (Es)c  for  an  arbitrary  S  G  G  ■  By  definition, 
Es  =  S(E).  Since  S  is  a  bijection  of  Rd  onto  itself, 

W  =  S(W  )  =  S(EU  Ec)  =  S(E)  U  S{EC) 


i.e.,  (S{E))C  =  S(EC).  So  we  have 

(Es)c  =  (S(E)Y  =  S(EC)  =  Si^Bj) 

Hence, 

H  =  (Hy  =  nseg  (EsY  =  ns&g  S^Bj) 
Notice,  if  5'  G  Q  and  L(S)  <  L(S'),  then 

S'OJJ^)  c  S'(u?=1^)  C  5(^5,). 

So, 
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H  =  nseg  *MU?) 
Therefore,  if  is  closed.  Since  it  is  also  bounded,  it  is  compact. 


Since  Q  =  U£L0£fc,  it  follows  from  Lemma  2.3.2  that 
Hc  =  US£g  Es  =  ur=0  U5eefc  Fs  = 


□ 


where  F0  =  F. 

For  any  5  €  5  -  {e},  if  S  =  InJn2  •  •  •  Ink-J*»  then  let  pi(5)  =  /nt,  the  first 
component  of  its  unique  representation.  For  example,  Pi(hhh)  =  h-  Hence,  for 
any  k  >  1, 

Fk  =  UsegkEs  =  U"=1  U5eet)p1(s)=/J-  F5 

where  n  is  the  total  number  of  inversions  defined  in  Rd  . 

Setting  F3k  =  USee,,F1(s)=/j^s  for  any  k  >  1,  we  have 

Ffc  =  U^Fi 

and 

tfc  =  ur=oFfc  =  F0  U  (U&  U?=1  Fjj). 
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This  partition  of  Hc  will  be  the  basis  for  defining  functions  on  Hc  in  the  next  chapter, 
and  ultimately,  the  reflected  diffusion. 

The  map  $  is  also  continuous,  as  shown  in  the  next  proposition.  The  proof 
of  this  proposition  is  similar  to  that  of  Proposition  2.11  in  [5],  except  that  since  the 
boundaries  of  Bj,  j  —  1,2,  ...,n  are  disjoint,  the  restriction  on  the  domain  of  $  in  [5] 
is  no  longer  necessary. 

Proposition  2.3.2.  $  is  continuous  on  Hc. 

Proof.  First  of  all,  by  the  definition  of  the  inversion,  any  member  of  Q  is  a  continuous 
map  on  Rd,  hence,  on  Hc. 

For  any  point  x  G  Hc,  either  there  exists  some  S  G  Q  such  that 

S(x)  G  U]=1dBj, 

or  there  is  no  such  S.  We  shall  discuss  these  two  cases  separately. 
First,  if  for  some  1  <  j  <  n  and  S  =  InJm  ■■■Ink, 

S{x)  G  dBj 

then 

x  G  S-'idBj) 

and 

=  $  o  S(x)  G  dBj 

since  the  boundaries  of  all  balls  are  disjoint  and  [x]  has  exactly  one  point  in  E. 
Because  5"1  o  $(x)  =  x,  there  exists  a  neighborhood  N  of  $(x)  such  that 

N  cEU{Bj-  Ij^jBi)). 

Since  $  is  the  identity  on  E  and  $  =  Ij  on  Bj  -  Ij(\Ji^jBi),  we  have 
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$\e(x)  =x 

and 

*l(fl>-/j(u^/B»))(:r)  =  cj  +  p^T^jp^  ~  c^ 
For  every  x  G  dBj,  \\x  -  Cj\\  —  1.  Then  we  have 

i.e.,  these  two  continuous  functions  coincide  on  dBj.  Therefore,  on  S~1(E),  $  =  5; 
on  S'1{Bj  -  Ij(UijijBi)),  $  =  IjS.  Both  of  these  functions  are  continuous  and  they 
coincide  on  $~l{dBj).  Therefore,  $  is  continuous  on  the  open  set  S~l(N)  containing 
x. 

Second,  if  for  any  S  G  Q , 

S(x)  $  \J0idBj, 

then  <J>(x)  £  U"=19Bj.  So  *(*)  €  £°,  the  interior  of  E.  Since  [a;]  has  exactly  one 
point  in  E,  there  exists  a  unique  T  e  Q  such  that  x  —  To  So         =  T_1(y) 

on  T(E°)  and  it  is  continuous. 

□ 


CHAPTER  3 

TIME  CHANGE  AND  //-TRANSFORM  FUNCTIONS 
3.1    The  Time  Change  Density  q 

Having  completed  the  partition  of  Rd,  we  are  now  ready  to  define  a  positive 
function  g  on  Hc.  This  function  g  will  then  define  a  time  change  of  the  Brownian 
motion. 

Recall  that 

H   =  USee*,Pi(s)=/;£s,    for  k  >  1,  1  <  j  <  n 
F0   =  E 

Fk   =   U^Fl  for  k  >  1 
Hc  =  U?=0Fk 

Notice  that  for  any  j  G  {1,2,...  ,  n}, 

Ij(Fl)   C  F0 

Ij(Fi)   c  Fk~i,  for  k  =  1,2,3,... 

In  the  following  definition,  we  shall  define  g  on  F0  first.  Then,  inductively,  we  can 
extend  the  definition  of  g  to  all  of  H°.  Let  us  now  state  the  definition  of  g  first  and 
then  prove  that  it  is  a  well-defined  positive  function.  Define 

g{x)  =  1,    x  £  F0; 

9(x)  =  ,,     1   „4g(/j(*)),  x  e  Ft-  j  €  {1,2, ...  ,  n},  *  =  1, 2,3, . . . 
IF  cj\\ 


(a) 
(b) 
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Proposition  3.1.1.  The  function  g  is  well-defined  and  positive. 

PROOF.  By  (a),  g  is  defined  on  F0.    For  all  x  e  Fx  =  U?=1F/,  if  x  e  F(, 

j  e  {1,2,3, . . .  ,n},  then 

/,(«)  e  f0 

and  by  (b), 

9{x)  -  w^¥9{Iiix)) 
i 

llz-cjll4 

So  g  is  well  defined  on  Fi  =  U"=1F/  and  it  is  positive. 

Suppose  that  g  is  defined  on  F0,  Fx,  Ffc_i  already.  Then  for  any  x  € 
F3k,Ij{x)  is  in  Fk-i,  on  which  g  is  already  defined.  So  g(x)  can  be  defined  according 
to  (b)  for  each  j  =  1,  2,  . . .  ,  n.  Therefore,  p  is  defined  on  Ffc  =  U^=1F^.  p  is  obviously 
positive  since  it  is  so  in  each  previous  inductive  step. 

□ 

The  role  of  g  will  become  clear  in  the  next  chapter,  where  it  introduces  a 
positive  additive  functional,  which  eventually  leads  to  a  time  change  of  the  diffusion 
process  we  study.   The  relationship  between  the  function  g  and  the  inversions  is 
demonstrated  in  Proposition  3.1.2  below. 
Proposition  3.1.2.  g(x)  satisfies 

for  each  j  =  1,2,  ...  ,n  and  all  x  G  Hc. 

Proof.  First,  for  every  x  6  F0  C  Hc,  Ij(x)  G  F{.  So,  g(x)  =  1  by  (a)  and 
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by  (b).  Therefore,  (3.1)  is  true. 

Second,  suppose  x  G  Fk  for  some  k  >  1,  then  x  G      for  some  1  <  i  <  n. 
There  are  two  cases. 

(i)  If  i  =  j,  then  Ij(x)  G  Fk-i-  According  to  (b),  the  left  hand  side  is  defined 
by  the  right  hand  side.  So,  the  above  identity  is  true. 

(ii)  If  i  ?  j,  then  Is(x)  G  Fi+V  By  (b)  , 


So,  the  right  hand  side 
Thus,  the  identity  is  proved. 


□ 


3.2    The  Function  h  and  Its  Properties 

We  now  define  a  superharmonic  function  h  which  will  be  used  to  /i-transform 
our  Brownian  motion. 

For  any  S  G  Q  -  {e},  let  as  =  Ds~l,  where  Ds  is  defined  as  in  Proposition 
2.2.2.  Also,  let 

C  =  {S{oo)  -.SeG}. 
Since  Q  is  countable,  so  is  C.  We  define 

Ead-2 
-  J.    ....        for  all  x  G  Rd  -  C. 
\\x  -  S(oo)||d-2 


seG  -{e} 
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Note  that 

lim  h(x)  =  oo 

x->z 

for  every  z  e  C,  so  we  may  set  h(z)  =  oo  for  every  z  e  C. 

It  is  obvious  from  this  definition  that  h  is  positive.  In  the  next  section,  we 
shall  prove  it  is  finite  in  Kd  except  on  the  countable  set  C. 

For  the  present,  let  us  examine  its  property  associated  with  inversions  first. 

Lemma  3.2.1.  For  S  ±  Ij,  Ds\\S{oo)  -  ij(oo)||  =  DIjS,  hence, 

_  as 
||5(oo)- 7,(00)11 

Proof.  There  are  two  possible  cases,  Pi(S)  =  Ij  and  Pi{S)  ^  Ij. 

(i)  If  Px (S)  =  Ij,  suppose  S  =  Ijln2  ...  Ink,  k>2.  Then  IjS  =  Jna  . . .  I„t,  so 


=   ||7,-7n2  . . .  7nfc(oo)  -  7,(00)11  x  y/^oo)  -  In2  . . .  /„t(oo)||  x 
||7„2(oo)-7n3  ...  /njk(oo)||  x--  -  x  p^oo)  -7„,(oo) || 
=   ||/„a(oo)-/„3  ...  7„fc(oo)||  x--  -  x         (00)  -/„, (00) || 
=  TJ/.s 

(ii)  If  Pi  (5)  ^  7„  suppose  S  =  In2  . . .  7nfc,  k  >  2.  Then  7,5  =  7,7n2  ...  7, 

||7,(oo)-S(oo)||Ds 
=   ||7,  (oo)  -  7„2  ...  7nfc(oo)||  x 

||/na(00)  -/„...  In*  (OO)  ||  X--  -  X  ||/nfc_1(oo)-/nt(oo)|| 
=  P>/.S 


□ 
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Proposition  3.2.2.  The  function  h  satisfies 

ftW=||x-/,|oc)||^^W)'  i  =  1,2, 

Proof.  By  the  definition  of  h, 

1  ^  af2 


n 


x-IjiooW-*  \\x-S(oo)\\«-2 
sey  -{ejj} 

1  ad~2 
h{Ij{x))  =  l+nrM_wM|,-2+  E 


||7,(x)- 7,(00)11^         ^  5(oo)||^ 

Recall  that  for  any  5  G  Q  -  {e},  Pi  (5)  is  the  first  inversion  in  the  unique  represen- 
tation of  5.  Using  (2.2)  for  the  second  term  in  the  above  expression  and  Proposition 
2.1.3  for  the  terms  in  the  summation,  we  have 

hMx))   =   l  +  \\x-  /jMH     +     J-  ||x_/.5(0o)||--»||5(oo)-/i(oo)||^ 

Se£  -{e./j} 

So, 

Hx-^-MII^^^^  ~  1+||x-/,(oc)||^+5e^e  ^       lie  -  7,-5(00)11^ 
Let 

£'  =  {S:  Se£  -{e,Ij},  f\(S)  =  /,-} 

and 

^*  =  {S:  Seg-{e,  /J,  Pi(5)^/,} 
Then  we  have  £  -  {e,  Ij}  =  GjuGj*.  Moreover, 

Qi  =  {I3S:  SeGj*} 

and 

GJ*  =  {IjS:  SeG3} 

By  Lemma  3.2.1, 

as 

OiljS 


||S(oo)-J,(oo 


So,  if  we  let  S'  —  IjS, 

L*  \\x  -7,5(00)11^ 

seG  -{e.Ij} 

-  _  /)S(oo)||'"  +  ^  ||i-/jS(oo)||"-2 

||x-5'(oo)||d-2  +  ^.\\x-S'(oo)\\d-2 
S'eG3'  S'eG3 


E 


d-2 


||x-5'(oo)|r-2 


Hence, 


-,  d-2 

1  1  as< 

\\x  -  I.iooW-^1^   =   1+\\x-  +  ||x-5'(oo)||^ 

by  definition. 

□ 

3.3    The  Superharmonicitv  of  the  Function  h 

It  remains  to  be  shown  that  the  function  h  we  constructed  in  the  previous 
section  is  indeed  superharmonic.  Notice  that  h  is  an  increasing  limit  of  superharmonic 
functions.  We  will  prove  that  this  limit  is  finite  except  on  a  polar  set. 
Proposition  3.3.1.  The  function  h  is  a  superharmonic  function. 
Proof.  Recall  the  construction  of  the  function  h: 

*(*)=i+  E  l|x-£)ir*-  **^~c- 


seG -{e} 
Now,  for  k  e  N,  let  us  define 


Se£  -{e},  L(S)<fe 
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Since  the  constant  function  and  functions  of  the  form 

/(a:)  =  p  -  fe||d"2 

are  all  harmonic  functions  in  Rd,  it  is  easy  to  see  that  for  any  k,  hk(x)  is  also  a 
harmonic  function.  Therefore,  it  is  superharmonic. 
Notice  that  hk+i(x)  >  hk(x)  for  all  k  and 

lim  hk(x)  =  h(x). 

fc->00 

Therefore,  h  is  an  increasing  limit  of  superharmonic  functions. 

By  the  unique  representation  of  S,  the  cardinality  of  the  set 

gk  =  {s\seg,  L(s)  =  k} 

is  n(n  -  l)fc_1  and  we  have 

G  -{e}  =  U?=1& 
Now  let  us  examine  each  individual  term  of  h. 

Since  for  any  i  ^  j,  ||q  -  Cj\\  >  n,  as  we  stipulated  in  chapter  2,  we  can 
estimate  the  magnitude  of  each  as- 

For  any  5,  Ds  is  the  product  of  L(S)  -  1  factors,  each  of  which  is  of  the  form 

117,(00)  -T(oo)|| 

where  T  6  Q  -  {e}  and  FX(T)  =  /<  ^  Ij.  Since 

T(oo)  e  Bi 

and 

/^(OO)  =  Cj 

it  follows  that 

||7,(oo)-T(oo)||>||Cl-cJ-l  =  n-l 
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Furthermore,  we  have 

inf{||/3-(oo)  -  T(oo)||,  j  =  1, 2, . . .  ,  n;  P1{T)  =  h  *  /,}  >  n  -  1 

or,  more  explicitly,  we  have 

5  =  -J—  inf {p^oo)  -  T(oo)||,  j  =  1,2, . . .  ,n;  i\(T)  -  U  ±  I,}  >  1 
n  —  1 

So,  for  any  S  such  that  L(5)  =  k,  Ds  >  {n  -  l)k-lSk~l.  Consequently, 

1 

as  <  ^(n-l)*-1 
For  any  given  x  e  E,  and  any  S  £  &  -  {e}, 

||x-5(oo)||  >  M 

for  some  M  >  0.  Thus, 


<f-2 


< 


and 


||x-5(oo)||d-2     (n_  i)(fc-i)(d-2)5(fc-i)(d-2)M<i-2 


^         Hx-5(oo)||d-2    <    (n  -  i)(fe-i)(rf-2)5(fe-i)(d-2)Md-2 


Se£  -{e},  L(S)=fc 


n  (    \  \k-l 


Md-2{n  -  l)(*-i)(rf-3)  V<^-2 

Hence,  the  infinite  series 


m*)  =  i+E  £ 


«f2 


fc=1  SeC? -{e},  L(5)=fc  " 

is  convergent  for  all  x  G  E,  i.e.,  /i(:r)  is  finite  in  E. 

Now  we  have  proved  that  function  h  is  a  superharmonic  function  which  is  finite 
except  on  a  polar  set. 

□ 


CHAPTER  4 
MAIN  RESULT 

4.1    The  Construction  of  X± 

In  Kd,  d  >  3,  the  inversion  of  Brownian  motion,  even  in  the  simplest  one-ball 
case,  is  no  longer  a  Brownian  motion,  as  is  demonstrated  in  [15],  p. 100.  It  is  an  h- 
transform  of  Brownian  motion  for  a  superharmonic  function  h.  The  question  is,  what 
if  we  start  with  a  process  which  is  already  an  /i-transform  of  Brownian  motion?  The 
answer  is  trivial  in  the  one-ball  case,  due  to  the  idempotent  property  of  inversions. 

Let  J  be  the  inversion  with  respect  to  the  unit  sphere,  centered  at  the  origin 
with  radius  1.  Suppose  Zt  =  I(Bt)  where  Bt  is  Brownian  motion.  Zt  then  is  an 
/i-transform  of  Brownian  motion,  where 

h{x)  =  w=~* 

It  is  easy  to  see  that  I{Zt)  =  I2{Bt)  =  Bt.  Notice  that  the  group  of  homeomorphisms 
on  Rd  introduced  by  the  inversion  is  finite,  namely,  Q  =  {e,  I}. 

However,  in  the  case  of  two  or  more  balls,  the  cardinality  of  this  group  be- 
comes infinite,  due  to  the  fact  that  inversions  with  respect  to  different  spheres  do  not 
commute.  This  means  that  the  idempotent  property  of  each  individual  inversion  will 
no  longer  give  us  a  clear  indication  as  to  what  the  Brownian  motion  will  "become" 
after  it  has  undergone  a  series  of  inversions. 

Our  solution  here  is  to  construct  a  diffusion  process  X  which  is  invariant  under 
any  of  these  inversions.  Specifically,  by  invariant  we  mean  that  if  5  is  a  member  of 
Q ,  then  for  any  x  G  Rd  -  C,  {Xt,  Px)  and  (S(Xt),  Ps-1(*))  are  identical  in  law. 

We  shall  show  that  this  particular  diffusion  can  be  obtained  by  an  /i-transform 
and  a  time  change  of  Brownian  motion. 
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Recall  that 

"(*)  =  '+  E  fi&f*'  forallI€Ri-c 

is  a  positive  superharmonic  function  and 

Rd  -C  =  {h<  00} 

Let  B  —  {Bt,T't)  be  the  standard  Brownian  motion  in  Rd  with  transition 
density  p(t,  x,  y)  and 

be  the  /^-transform  of  B.  The  transition  density  of  Bh  is  then 

nh(f  -r  ,a  _,/  h(y)p(t1x1y)h(x)-1   if  (x.y)  6  (Rd  -  C)  x  Rd 

p  j  o  if  (x)2/)eCxRd 

Then  P71  is  a  Markov  process  with  state  space  Rd,  transition  density  ph(t,x,y),  and 
initial  distribution  supported  by  Rd  —  C.  Bh  also  has  an  almost  surely  continuous 
modification  and  the  strong  Markov  property. 

Let  Ct,(Rd)  be  the  class  of  bounded  continuous  functions  on  Rd.  If  Pt  is  the 
semigroup  of  B  and  P/1  is  the  semigroup  of  Bh,  then  for  any  /  €  C6(Rd),  £  >  0  and 

Pthf(x)  =  -^-)Pt(hf)(x)  (4.1) 

In  Rd,  d  >  3,  it  is  easy  to  see  that  C  is  a  polar  set  of  B  since  it  is  countable. 
By  Doob  [12],  p.674,  almost  no  path  of  Bh  hits  C  at  a  strictly  positive  time. 
Next,  we  consider  a  time  change  of  Bh. 
Define 


At=  f  g{Bhs)ds, 
Jo 

where  g  is  the  positive  continuous  function  defined  in  section  3.1.  Hence  At  defines 
a  perfect  continuous  homogeneous  additive  functional  (PCHAF)  as  in  [18].  Let  r(t) 
be  its  right-continuous  inverse,  i.e., 

r(<)  =  inf{s  :  A,  >  t] 


so 

{r{s)  <t}  =  {At  >  s} 

Since  {J7?}  is  right-continuous,  for  each  s,  r(s)  is  an  {T^}  stopping  time.  Notice 
that  because 

g(x)  >  i 

for  all  x  G  H°,  we  have 

As  >  s 

for  all  s  >  0.  Hence,  for  all  t  >  0, 

r(t)  <  t  (4.2) 

This  is  a  useful  condition  in  formulating  the  martingale  problem  later. 

Set  Xt  -  B*(t)  and  Tt  =  T*{t).  Then  X  -  {XuTt)  is  a  time  change  of  Bh. 
This  means  that  X  is  a  strong  Markov  process  just  as  Bh  is.  The  time  change  of  Bh 
preserves  the  continuity  of  its  paths;  so  X  also  has  continuous  paths. 

In  conclusion,  starting  from  a  Brownian  motion  in  Rd,  we  have  constructed 
a  diffusion  process  X  path-by-path  in  Rd  —  C.  In  the  next  three  sections,  we  shall 
characterize  its  distribution  and  then  examine  its  invariance  and  symmetry  properties. 

4.2    Infinitesimal  Generators 

Having  constructed  the  diffusion  process  X  path-by-path,  we  can  now  examine 
its  distribution.  We  can  achieve  this  by  characterizing  its  infinitesimal  generator. 

Let  G  be  the  infinitesimal  generator  of  Brownian  motion  in  Rd  and  Vq  be  its 
domain.  Let  Co(Rd)  be  the  space  of  bounded  continuous  functions  on  Rd  which  vanish 
at  infinity.  Let  A  be  the  Laplace  operator.  The  natural  domain  of  the  operator  |A 
is 

=  {/  €  C0(Rd)  :  l-M  exists  and  ±Af  G  C0(Rd)} 
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For  d  >  2,  the  operator  G  is  the  closure  of  |A  in  the  following  sense:  /  6  £>G  if 
and  only  if  there  exist  functions  /„  in  T>(^A)  and  a  function  g  in  C0(Rd)  such  that 
||/„  -  /||  ->  0  and  ||±A/„  -  g\\  ->  0;  and  then  Gf  =  g. 

Let  Co(Rd  —  C)  denote  the  class  of  twice  continuously  differentiable  functions 
with  compact  support  on  Rd  —  C.  Then  clearly 

C20(Rd  -C)c  V(±A) 

By  the  definition  of  h,  it  is  twice  continuously  differentiable  on  Rd  —  C.  So, 
for  any  /  G  C$(Rd  -  C),  we  have 

hf  6  C02(Rd  -  C)  C  2>(^A) 

From  this  we  can  compute  the  infinitesimal  generator  of  P*,  the  semigroup  of 
Bf.  Recall  equation  (4.1)  and  let  G'1  denote  the  infinitesimal  generator  of  P/1.  Then 

Pthf(x)  -  f{x) 


Ghf{x)   =  lim 


t44-0  t 
1   ,im  Pt(hf)(x)  -  h(x)f(x) 


h(x)  mo  t 

It  is  necessary  to  point  out  here  that  in  Rd  —  C,  d  >  3,  the  domain  of  Gh  is 
richer  than  Cl{Rd  —  C).  What  has  been  established  here  is  that  P/1,  the  semigroup 
of  Bf,  has  infinitesimal  generator  Gh  where 

Ghf  =  ^A(hf) 

for  /  €  C02(Rd  -  C). 

Next  we  turn  to  investigate  the  infinitesimal  generator  of  the  semigroup  of  Xt 
by  studying  the  related  martingale  problem. 

From  the  above  discussion,  we  know  that  for  any  /  G  Co(Rd  —  C), 


C{  =  f{Bl)-  fGhf{Bhs)ds 
Jo 
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is  a  Px-martingale  for  any  x  G  Rd  -  C.  By  the  definition  of  r(t)  and  equation  (4.2), 
we  know  that  r(t)  is  an  optional  stopping  time  with  respect  to  JF/1  and 

Cl(t)  =  /Wfc))  ~  /  Ghf(Bi)ds 

J  0 

is  also  a  P^-martingale  for  any  x  G  Kd  —  C  with  respect  to  T^. 
Letting  s  =  t(u),  we  have 

u  =  As=  f  g(Bhv)dv 
Jo 

It  then  follows  that  du  =  g(B^)ds  and 


=  /(*•)-  /  g-'(X,)Gkf(Xu)du 

JO 

Since  r(0)  =  0, 

ExCfT{t)  =  ExCfT{0)  =  E*Cl  =  E*f(Bh0)  =  f(x) 

and 

f{x)   =  E*f(Xt)-  f '  E*{g-\Xu)Ghf{Xu))d'< 

Jo 

=  Ptf(x)-  f  Pu{g-lGhf){x)du 
Jo 

where  Pt  is  the  semigroup  of  Xt. 
Now  we  have 

Ptf(x)  -  f(x)  =  f  Pu{g-lGhf){x)du 

and  if  we  denote  the  infinitesimal  generator  of  Pt  by  >4, 

Ptf(x)  -  f{x) 


A  fix)   =  lim 


tuo  t 

{g-lGhf)ix) 

1  :Ghf(x) 


-—\—-Aihf)ix) 
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Thus,  we  have  established  that  the  infinitesimal  generator  of  Pt,  the  semigroup 
of  X,  acting  on  C02(Rd  -  C)  is 

Again,  C\  (Rd  —  C)  is  only  a  subspace  of  the  domain  of  A.  The  latter  is  in  fact 

richer. 

4.3    Invariance  of  Xt  and  The  Horizon  H 

In  this  section,  we  shall  show  that  the  process  Xt  is  invariant  in  law  under 
inversion.  Based  on  this,  we  can  further  prove  that  the  horizon  set  H  is  a  polar  set 
for  Xt. 

Theorem  4.3.1.  For  any  5  E  G ,  the  process  (S(Xt),  Ps~1(-x))  and  the  process 
(Xt,  Px)  are  identical  in  law. 

Proof.  Our  proof  is  based  on  the  fact  that  if  the  potential  of  (S(Xt),  P3'1^)  is  the 
same  as  that  of  (Xt,  Px),  then  they  are  identical  in  law. 

It  is  sufficient  to  prove  the  case  of  S  =  Ij  for  any  1  <  j  <  n.  In  the  following 
proof,  let  Ex  denote  the  expectation  with  respect  to  the  distribution  of  Xt  starting 
at  point  x  and  let      denote  that  of      starting  at  point  x. 


where 


E 
E 


poo 

Em  /  f(Ij(Xt))dt 
Jo 

poo 

tM  /  f&(BhTt))dt 
Jo 

poo 

/,(x)  /  f(Ii{Bt))g{Bi)da 
Jo 

u}(  vx  [.  I  Pt(Ij(x),  y)f{Ij(y))g{y)h(y)dsdy 
h{Ij{x))  Jw1  Jo 
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Through  an  elementary  calculation,  if  we  denote  |r(^  —  l)ir  2  by  K,  we  have 


poo 

/  Pt(I](x),y)dt  =  K\\y-Ij(x)f 
Jo 


Therefore,  if  we  let  z  —  Ij(y),  y  =  Ij(z),  and  recall  that  lP|  =  \\z  —  Cj\\  2d, 
we  have 


poo 

EI]{X)  /  f{Ii{Xt))dt 
Jo 


Recall  that 

g(Ij(z))  =  \\z-cj\\*g(z) 

by  Proposition  3.1.2  and 

h(Ij(z))  =  \\z-cj\\d-2h(z) 
by  Proposition  3.2.2.  Furthermore,  we  have 

Ik -2/11 


\\Ii(x)  -  Ij(y)\\  = 


according  to  Proposition  2.1.2.  Therefore, 


poo 

Bm  \  f(Ij(Xt))dt 
Jo 


K 

x 


h(x)\\x  -  Cj\\*-> 

JRd  f{z)g{z)\\z  -  Cj\\%{z)\\z  -  i^J^^jpj*  -  c3\\-2ddz 
After  simplifying  the  right  hand  side,  we  get 

pOO  TS  p 

EIj(x)  J  nij(Xt))dt=—J^df(z)g(z)h(z)\\z-x\\2-ddz 


43 


On  the  other  hand, 

poo 

Ex  f(Xt)dt 
Jo 


poo 

Ehx  f(BhTt)dt 
Jo 

poo 

=  Ehx  f{Bhs)g{Bhs)ds 
Jo 

=  W)LL  ^^^y)h{y)dsdy 

=  Ldf(y)9(y)h(y)\\y-x\\2-ddy 

Thus,  these  two  processes  are  identical  in  law  by  Blumenthal  and  Getoor  [7]. 

□ 

Now  consider  the  horizon  set  H.  As  shown  in  section  2.3, 

H  =  nseg  SiUj^B,) 

For  any  S  G  Q  —  {e},  let  t(S)  be  the  index  of  the  rightmost  inversion  in  the 
unique  representation  of  S,  i.e.,  for 


9  —  T   T  T 


let  t(S)  =  nk. 

Then,  for  any  n  >  0,  let 

where  Gn  =  {S  :  S  e  G ,  L(S)  =  n}.  It  is  clear  that  H  C  Hn+1  C  Hn  for  all  n. 
Let 

T  =  inf{* 

and 


Tn  =  inf{*  :  Xt  e  Hn}, 
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then,  since  i/n+i  C  Hn, 

{Tn+1  <  00}  C  {Tn  <  00} 

and 

{T  <  00}  =  n~  x{rB  <  00} 

We  need  the  following  lemma  to  prove  that  H  is  a  polar  set  for  Xt. 
Lemma  4.3.1.  For  any  x  e  Uj=1dBj, 

PX(T2  <  00)  <  m  <  1 

Proof.  First,  since  Xt  is  transient  and  H2  is  compact, 

PX(T2  <  00)  <  1 

for  any  x  G  Wj=ldBj. 

Since  x  h>  Px(r2  <  00)  is  harmonic  on  the  interior  of  H2,  it  is  continuous  on 
the  compact  set  U"=13Pj.  Therefore,  by  a  standard  argument  in  calculus, 

m  -     max    PX(T2  <  00)  <  1. 
xeu*mldBj 

□ 

Now  we  are  ready  to  prove  the  following  theorem. 
Theorem  4.3.2.  For  any  x  G  Rd  -  C, 

PX(T  <  00)  =  0 

Proof.  For  any  n  >  1,  x  e  Rd  -  C, 

Px(Tn+1  <  00)   =   Px(Tn  +  THn+l  o  6Tn  <  00) 

=  pa:(rn<oo,p^(r/fn+1  <oo)) 

=   Px(PXTHTHn+1  <  °o)Px(Tn  <  00)) 


For  every  u  €  J2,  Xr„(w)  G  S(U]=1dBj)  for  some  5  €  £„_i.   Hence,  by 
Theorem  4.3.1, 

PXTn(THn+l  <  oo)   =  P*T«(*t  €  //„+i  for  some  £  <  oo) 

=  Ps~^XT«\S{Xt)  e  Hn+l  for  some  t  <  oo) 

=  Ps~*<XT»>(J»f,  €  S_1(ffn+i)  for  some  *  <  oo) 

<  P5_1(x^)(Xt  e  H2  for  some  t  <  oo) 

Since  S_1(Xt„)  6  Vnj=ldBj,  by  Lemma  4.3.1, 

PXT"{THn+1  <oo)  <ro<  1. 

Thus, 

P*(Tn+a  <  oo)   <  mPx(Px(Tn<<x>)) 
=   mPx{Tn  <  oo) 

By  induction  on  n,  we  have 

Px(Tn+1  <oo)<mn  for  all  n  >  1. 


Therefore, 


P*(T  <  oo)   =    lim  Px(Tn+1  <  oo) 


<    lim  mn 

n—too 


for  all  x  e  Rd  -  C. 
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A  set  theory  analysis  shows  that  \H\  =  2  for  the  case  n  =  2.  For  n  >  2,  H  is 
a  Cantor-like  set. 


For  n  =  2, 


where 


and 


H  =  H1U  H2, 


H^n^ihhrhiBi), 


H2  =  n£o(/2/i)B/2(Bi). 

Both  i/i  and  i/2  are  intersections  of  nested  closed  sets  in  Rd.  Since  the  radii  of 
(IiI2)nIi(B2)  and  (hh)nh{Bi)  approach  zero  as  n  -»  00,  each  of  i?i  and  #2  contains 
exactly  one  point,  respectively.  So,       =  2. 

For  n  >  3,  rewrite  H  as 

I  =  h*!  uL(s)=n  5(u^f(s)^) 

Then  we  can  see  that  i/  is  the  intersection  of  a  sequence  of  nested  closed  sets.  There- 
fore H  is  closed. 

4.4    Symmetry  of  Xt 

In  this  section,  we  briefly  summarize  some  important  facts  about  dual  pro- 
cesses and  symmetric  processes  relevant  to  our  study.  It  then  follows  that  the  strong 
Markov  process  Xt  is  a  symmetric  process. 

Suppose  (fi,  Zt,  Pt)  is  a  Markov  process  in  Rd  and  U  is  its  potential  operator. 
Then,  for  any  /  <E  C°°(Rd),  we  have 


Uf(x)=  r  Ptf{x) 
Jo 


dt. 
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Suppose  the  semigroup  operator  Pt  has  a  density  function  Pt(x,  y)  with  respect 
to  a  measure  £,  then  it  is  easy  to  see  that  U  has  a  kernel  u(x,  y)  with  respect  to  £ 
and 

Uf(x}  =  Jf(y)u(x,y)ady), 

where 

poo 

u(x,y)  =  /  pt(x,y)dt. 
Jo 

The  dual  process  of  (f2,  Zt,  Pt)  with  respect  to  77  is  defined  to  be  a  Markov 
process  (Q,  Zt,  Pt)  such  that 

J^d  Ptf(x)g(x)r](dx)  =  f(x)Ptg(x)V(dx) 

If  U  is  the  potential  operator  of  (Q,,Zt,Pt)  with  kernel  u(x,y)  —  u(y,x),  then 
we  have 

[  (Uf(x))g(x)r,(dx)  =  [  f(x)U(g(x)Hdx). 

If  U  =  U,  we  say  Markov  process  (fi,  Zf ,  Pt)  is  symmetric  with  respect  to  77. 

Next,  let  us  consider  the  /i-transform  of  (fl,Zt,Pt),  for  some  positive  super- 
harmonic  function  h.  If  we  denote  this  process  as  (fi,  Z£,  P/1)  and  let  (£l,Zt,Pt)  be 
the  dual  process  of  (fi,  Zt,  Pt)  with  respect  to  reference  measure  C,  then  we  have  the 
following  fact. 

Proposition  4.4.1.  Suppose  (Q,  Zt,Pt)  is  a  Markov  process  and  Pt  has  density 
function  pt {x,  y)  with  respect  to  (  and  (Q,  Zt,  Pt)  is  its  dual  process  with  respect  to 
a  reference  measure  C-  Then  for  any  positive  superharmonic  function  h,  the  Markov 
processes  (Q,  Zf,  P*)  and  (Q,  Zu  Pt)  are  dual  processes  with  respect  to  the  reference 
measure 

Proof.  We  know 
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So,  for  any  /,  g  €  C°°(Rd)  , 

L(P?m)g(x)h(xK(dx) 
=  JRd  [JRd  /Gflrt  (*,  y)ady)]g(x)h(x)((dx) 

=   JRd  [JRd  f(y)h(y)Pt(x,  y)((dy)}g(x)((dx) 

=  JRd  f(y)h(y)iJRd  Mx>  y)g(x)((dx)}ady) 

=   JRd  f(y)iJRd  Mv,  x)g(x)C(dx))h(y)ady) 

=  JRdf(y)(Pt9(y))HyK(dy) 

Hence,  (Cl,  Z?,Pf)  and  (f),  ZuPt)  are  dual  to  each  other  with  respect  to  hC,. 

□ 

Now  that  we  have  established  that  (0,,  Zf,  Pf)  and  (ft,  Zt,  Pt)  are  dual  pro- 
cesses, we  may  use  this  result  once  more  to  conclude  the  following  result. 
Corollary  4.4.1.  The  Markov  processes  (f2,  Z£,  P/1)  and  (ft,  Zf,  Pf)  are  dual  pro- 
cesses with  respect  to  h2C,. 

Proof.  By  Proposition  4.4.1,  (ft,  Zt,Pt)  is  the  dual  process  of  (ft,  Z^P/1)  with  re- 
spect to  h(.  So,  since  the  /i-transform  of  the  former  is  (ft,  Z*,  Pth),  using  Proposition 
4.4.1  again,  we  conclude  that  this  process  is  the  dual  process  of  (ft,  Z^Pf1)  with 
respect  to  h2(. 

□ 

Now  we  are  ready  to  state  the  central  result  of  this  section. 

Theorem  4.4.1.  If  a  Markov  process  (ft,  Zt,Pt)  is  symmetric,  then  so  is  its  h- 
transform  (ft,  Z*,  P/1). 
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Proof.  Suppose  (ti,  Zt,Pt)  is  symmetric  with  respect  to  (.  Then  (Q,Zt,Pt)  = 
(n,ZuPt).  So,  (ft, Pf)  =  (Q,Z{>,Pth).  By  definition,  {Q,Z£,Pf)  is  symmetric 
with  respect  to  h?(. 

□ 

From  this,  we  have  the  following. 
Proposition  4.4.2.  The  Markov  process  Xt  is  symmetric. 

Proof.  In  fact,  from  the  discussion  above,  we  can  identify  the  potential  V  of  (ft,  B*,  P£) 
for  any  /  G  C°°(Rd)  and  x  G  Hc, 

Vf(x)  =  J   f{y)v(x,  y)v(dy) 

where  rj  —  h2X  and 

f00  p(t,  y) 

70  h{x)h(y) 

Let  £/  denote  the  potential  of  Xt.  Let  i?x  and       denote  expectations  with  respect 
to  the  distribution  of  Xt  and  that  of  Bf,  respectively.  Then  for  any  /  G  C00^) 
and  x  G  ifc, 

/•oo 

Uf(x)   =   Ex  /  /(Xt)dt 
Jo 

poo 

=  Ehx  }{B\t))dt 
Jo 

poo 

=  Ehx  f(Bi)g(Bt)d8 
Jo 

=  V(gf)(x) 

=  J^df(y)v{x,y)g{y)v(dy) 
=  J^df(y)v(x,y)Z(dy) 

where  £  =  grj  and  v(x,y)  is  symmetric  as  shown  above.  Hence,  Xt  is  symmetric. 

□ 
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By  now,  we  have  proved  that  Xt  is  a  symmetric  diffusion  process  which  is 
invariant  in  law  under  any  S  £  Q . 

4.5   The  Reflected  Diffusion  Wt 

Now  we  are  ready  to  identify  the  reflected  diffusion  process  in  E.  Recall  that 
$  is  continuous  from  Rd  onto  E.  It  has  the  property 

for  every  S  E  Q .  This  means  that  $  is  constant  on  any  orbit  [x]. 

Lemma  4.5.1.  For  any  continuous  function  /  on  E,  Pt(f  o  <fr)(x)  is  constant  on  the 

orbit  [x]. 

Proof.  By  Theorem  4.3.1,  the  processes  {Xt,Px)  and  (S(Xt),  Ps_I(l))  are  identical 
in  law  for  any  S  G  Q .  Thus,  for  any  continuous  function  /  on  E, 


Px[(fo^)(Xt)}  =  Ps-lx[(fo^(S(Xt))) 
But  since  $o,S  =  $,/o$oS,  =  /o$.  We  have 

Px[(fo^)(Xt)]  =  Ps-lx[(fo^)(Xt)} 

Therefore, 


Pt(fo<S>)(x)  =  Pt(fo<S>)(S-1(x)) 
for  any  S  e  Q .  This  proves  that  Pt(/  o  $)(a;)  is  constant  on  the  orbit  [x]. 


□ 


Corollary  4.5.2.  Let  {Ua;a  >  0}  be  the  resolvent  of  the  semigroup  {Pt;t  >  0}. 
Then  for  every  continuous  function  /  on  E, 


poo 

Ua(f  o         =  Ex  /    e-Qt(f  o  *){Xt)dt 
Jo 


is  constant  on  [x]. 

Thanks  to  these  facts,  we  can  define  a  Markov  transition  function  on  (E,  13(E)). 
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Definition  4.5.3.  For  any  t  G  [0,  oo),  x  G  E  and  A  G  B(E),  let 

Qt(x,A)  =  P?(Xte$-1(A)) 

where  ye  [x]. 

By  Lemma  4.5.1,  Qt{x,  A)  is  well  defined  since  Pt(f  o  $){x)  is  constant  on  the 
orbit  [x].  The  next  theorem  shows  that  it  is  a  Markov  transition  function. 
Theorem  4.5.4.  The  collection  {Qt(-,  -),t  6  [0,  oo)}  is  a  Markov  transition  function 
on  (E,  B(E)). 

Proof.  First,  since  Xt  is  a  Markov  process,  for  any  x  G  E  and  A  e  B(E), 
A  ->•  Qt(x,A)  is  a  probability  measure  and  x  ->•  Qt{x,  A)  is  B(£)-measurable. 
Second,  if  we  write 


Qtf(x)=  I  f(z)Qt(x,dz) 
Je 


then 


Qt(x,A)  =  (QtIA)(x) 

=  (Pth-HA))(y) 
=  Pt(lAO$)(v) 

for  any  y  €  [x]. 

By  a  standard  measure  theory  argument,  for  any  bounded  continuous  function 
/  on  E,  we  have 

Qtf(x)  =  Pt(fo  *){y) 

for  any  y  G  [x]. 

Therefore,  for  any  s  >  0, 

Qt+s/(*)  =        °  *)(y) 

=  Pt(Fs(/o$))(y) 
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Since  Ps(f  o  $)(z)  =  Qsf{x)  for  x  =  $(z),  we  have 

Ps(/o$)(z)  =  ((Qs/)o$)(z) 

therefore, 

Qt+Sf(x)   =  Pt((QJ)o$)(y) 
=  Qt(Qsf)(x) 

Thus,  {Qt(-,  -))^  ^  [0,oo)}  is  a  Markov  transition  function. 

□ 

We  can  also  define  a  resolvent  on  E  by  setting 

Vaf(x)  =  Ua(fo$)(y) 

where  y  €  [x]. 

For  every  u>  in  the  underlying  probability  space  (fi,  T,  P),  define 

Wt(w)  =  $(Xt(w))  for  all  <  >  0. 

Notice  here  that  a  path  of  Wt  is  derived  from  the  path  of  Xt.  It  is  easy  to  see  by  the 
definition  of  $  that  for  every  fixed  u  the  path  <&(Xt(u)))  lies  in  E. 

The  next  theorem  shows  that  Wt  is  a  Markov  process  in  E  with  Qt  as  its 
transition  function. 

Theorem  4.5.5.  Wt  is  a  strong  Markov  process  in  E  with  continuous  paths. 
Proof.  First,  by  the  above  construction,  all  paths  of  Wt  lie  in  E  and  they  are  all 
continuous. 

For  every  bounded  continuous  function  /  on  E  and  every  x  G  E, 

Qtf(x)    =  Pt(/o*)(ar) 
=  E*(/o$)(Xt) 
=  Exf(Wt) 
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where  Ex(-)  is  the  expectation  with  respect  to  the  distribution  of  Wt  starting  at  x. 

Moreover,  for  any  bounded  continuous  function  /,  the  map  t  — >  Vaf(Wt)  is 
right  continuous  since 

Vaf(Wt)  =  UQ(fo$)(Xt). 

By  the  well  known  theorem  of  Blumenthal  and  Getoor,  (  [7],  Chapter  1,  Theorem 
(8.11)  ),  Wt  is  a  strong  Markov  process. 

□ 

By  the  definition  of  Qt,  the  Markov  process  Wt  also  inherits  other  properties 

oiXt. 

Theorem  4.5.6.  Wt  is  a  symmetric  process. 

Proof.  If  Pt  is  the  dual  semigroup  of  Pt,  then  Pt(f  o  $)(x)  is  also  constant  on  [x]. 
So  we  can  define  Qt  similarly,  and  it  is  the  dual  semigroup  of  Qt.  Since  Pt  =  Pt  by 
the  symmetry  of  Xt,  we  have 

Qt  =  Qt- 

Therefore,  Wt  is  a  symmetric  process. 

□ 

Theorem  4.5.7.  Wt  is  a  transient  Markov  process  with  infinite  lifetime. 
Proof.  Let  K  C  E  be  a  compact  set  such  that 

\{K)  >  0 

where  A  is  the  Lebesgue  measure.  Then  K  is  bounded  and  closed.  By  the  definition 
of  is  also  bounded  and  closed.  Therefore,  is  compact  in  Rd  and 

\($-1{K)>0). 

By  the  transience  property  of  Xt,  for  any  x  E  E, 
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o  <  Uh-i(K)(y)  <  oo 

for  every  y  £  [x].  Since 

VIK(x)  =  UIKo*(y)  =  Uh-HK)(y), 

it  then  follows  that 

0  <  VIK(x)  <  oo. 

Thus,  Wt  is  also  transient. 

The  fact  that  h  is  a  polar  set  for  Xt  implies  that  Xt  has  infinite  lifetime  in  Hc. 
Hence,  the  lifetime  of  Wt  in  E  is  also  infinite.  □ 

Therefore,  we  have  demonstrated  that  Wt  is  a  symmetric,  transient  diffusion 
process  with  infinite  lifetime  on  E. 

Theorem  4.5.8.  Wt  is  an  /i-Brownian  motion  in  the  interior  of  E  and  is  normally 
reflected  on  the  boundary 

dE  =  Unj=xdBj 

Proof.  Let  Cq(E)  be  the  class  of  twice  continuously  differentiable  functions  with 
compact  support  on  E.  Let  B  be  the  infinitesimal  generator  of  Qt.  Then,  for  every 
/  6  C$(E),  x€E, 

JK  '     tuo  t 

and 

f(Wt)  -  f{Bf){Wa)ds 
Jo 

is  a  P^-martingale  for  every  x  E  E.  Since 

Wt(u)  =  $(Xt(u))  for  all  t  >  0. 

we  know  that 


/o$(Xt)-  [\Bf)o^(Xs 
Jo 


)ds 


55 

is  a  P^-martingale.  Hence,  we  have 

Pt(f  o  $)(«)  -  /  o  *(*)  =  r  P4((B/)  o  $)(x)ds 

Therefore, 

(Bf)o*(x)  =  Hmi/  £((B/)o 

1440  I  Jo 

=   M±(£(/o*)(s) -/<,*(*)) 
=  4(jfo#)(x) 

For  any  a;  G  there  exists  a  neighborhood  of  a;,  jV(x),  such  that  N(x)  C  i?". 
Therefore,  we  have 

$(z)  =  z 

for  all  z  G  N(x).  Thus,  since  g(z)  =  1  for  z  E  E,  it  follows  that 

Bf(x)  =  Af(x)  =  ^&(hf)(x) 

This  shows  that  Wt  is  an  /i-Brownian  motion  in  E°. 

If  x  G       =  Uj=1dBj,  then  there  exists  a  neighborhood  of  x,  7V(a;),  such  that 
N{x)  C  BjCiE  for  some  j.  For  any  z  G  iV(x)  n  £, 

$(z)  =  * 

and  for  any  z  G  A^(x)  D  Bj, 

k*)  *  %{*). 

Let  the  center  of  Bj  be  Cj  —  (cpdj, . . .  dj).  Let  f(z)  =  f(z\,  z2,...  ,  £</)  and 

_  .df_  df_  df_ 
dz\ '  dz2 '  ' 

Compute  V(h(f  o  J^)),  we  get 

v(M/ <>/,■))(*)  =  (Vft)(«)(/o /,-)(«) +  M«)v(/o /,•)(«) 


=  (VAK*)/(W)+*(*)(V/K/iW) 


|z-c,||2  llz-c-114" 
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where  J  is  the  d  x  d  identity  matrix  and 


D(z,  Cj) 


z2-c 


(zi  -  c],z2-  c2,...  ,zd-cj) 


Compare  this  with 


V(hf)(z)  =  Vh(z)f(z)  +  h(z)Vf(z) 


Notice  that  x  G  dBj  implies  Ij(x)  =  x,  \\x  —  Cj\\  =  1,  and  with  g(x)  =  1,  we  have 


Bf{x)  = 


h(x) 


(Vf)(x)D(x,Cj) 


for  every  x  e  dBj,  j  —  1, 2, . . . ,  n. 

Notice  that  for  every  j  and  any  vector  v  =  (v\,  v2, . . .  ,  Wrf), 

vD(x,  Cj)  =  K(xi  -  c),  x2  -  c2j, . . .  ,xd-  cj) 

for  some  scalar  K.  So  vD(x,Cj)  is  normal  to  the  surface  dBj. 

From  the  symmetry  property  of  Wt  and  the  fact  that  B  is  normal  to  dE 
on  every  boundary  point,  we  know  that  Wt  is  an  /i-Brownian  motion  with  normal 
reflection  on  the  boundary. 

□ 


This  concludes  the  construction. 


CHAPTER  5 
CONCLUSION 

5.1  Summary 

In  the  study  of  multidimensional  diffusion  processes,  there  are  generally  two 
approaches,  as  discussed  in  Stroock  and  Varadhan  [8].  The  analytic  approach  is  one  of 
them.  In  this  approach,  a  diffusion  process  is  studied  through  the  parabolic  equation 
associated  with  it,  such  as  in  the  study  by  Fukushima  [9].  The  major  difficulty 
with  this  approach  is  that  it  is  generally  a  non-trivial  problem  to  prove  that  this 
parabolic  equation  has  a  sufficiently  large  family  of  solutions  with  required  properties 
to  characterize  the  associated  diffusion  process.  And  since  we  can  only  learn  about 
this  process  by  studying  the  associated  parabolic  equation,  probabilistic  intuitions 
are  somewhat  obscured  by  analytic  techniques. 

Another  approach  is  to  study  diffusion  processes  using  Ito's  stochastic  inte- 
gration. Examples  along  this  line  can  be  found  in  Ikeda  and  Watanabe  [10]  and 
Watanabe  [11].  Since  the  diffusion  processes  are  studied  in  association  with  the  un- 
derlying Brownian  motion,  more  probabilistic  insight  can  be  gained.  In  particular,  a 
process  can  be  constructed  path- by-path,  where  the  state  space  of  the  process  and 
the  underlying  probability  space  are  explicit. 

Glover  and  Rao  gave  a  path-by-path  construction  of  reflected  Brownian  motion 
in  R2  (See  Glover  and  Rao  [5]).  They  achieved  this  by  geometrically  reflecting  the 
paths  of  a  standard  Brownian  motion  on  the  boundary  of  an  open  set  and  showing 
that  the  resulting  process  is  a  Brownian  motion  with  a  time  change. 

In  this  thesis,  a  reflected  h-  Brownian  motion  in  Rd,  d  >  3  is  constructed  path- 
by-path.  The  definition,  properties  of  /^-Brownian  motion  and  its  association  with 
Brownian  motion  are  discussed  in  Doob  [12]. 
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This  path-by-path  construction  starts  with  a  standard  Brownian  motion  in 
Rd,  d  >  3.  An  /i-Brownian  motion  is  then  derived  from  it  for  a  superharmonic  function 
h.  This  /i-Brownian  motion  is  then  time-changed  to  a  diffusion  process  Xt.  Under 
the  conditions  imposed  on  the  superharmonic  function  h  and  the  time  change,  Xt  is 
a  symmetric  Markov  process  with  infinite  lifetime.  A  continuous  many-to-one  map 
$  :  Rd  ->  E  is  defined,  where  E  =  Rd  -  Unj=lBj. 

It  is  then  proved  that  Wt  =  ${Xt)  is  a  strong  Markov  process  with  state  space 
E.  Wt  is  also  symmetric  and  has  infinite  lifetime.  The  infinitesimal  generator  of  this 
process  is  then  computed  and  it  shows  that  Wt  behaves  like  /i-Brownian  motion  in 
the  interior  of  E  and  has  normal  reflection  on  the  boundary. 

It  is  worth  noticing  that  a  many-to-one  function  in  general  does  not  map  a 
Markov  process  to  a  Markov  process.  A  study  in  this  subject  can  be  found  in  Glover 
and  Mitro  [6]. 

5.2  Discussion 

The  interest  in  studying  reflected  /i-Brownian  motion  originated  in  the  study  of 
reflected  Brownian  motion  and  Markov  processes  in  general.  The  connections  between 
Markov  processes  and  martingales,  probability  theory  and  potential  theory,  provide 
powerful  tools  in  the  study  of  these  subjects.  References  on  these  subjects  can  be 
found  in  Chung  [13],  Rao  [14],  Durrett  [15],  Freedman  [16],  Revuz  [17],  Williams  [18], 
McKean  [19],  Ito  and  McKean  [20]  and  Landkof  [21]. 

As  pointed  out  in  Chung  [13],  an  analytical  approach  to  the  reflecting  boundary 
problem  in  multidimensional  spaces  presents  a  major  difficulty.  In  the  case  studied 
here,  the  construction  is  realized  path-by-path. 

There  are  three  restrictions  imposed  on  this  construction,  namely,  the  number 
of  open  balls  Bj,  their  radii,  and  the  minimum  distance  between  them. 

The  restriction  that  each  ball  has  radius  1  is  for  the  sake  of  clarity  only  and 
can  be  removed  without  altering  the  the  basic  argument.  That  is,  this  construction 


will  not  change  fundamentally  if  the  balls  are  allowed  to  have  different  radii,  but  the 
symbols  and  formulas  will  be  more  complicated. 

The  number  of  balls  and  the  minimum  distance  between  them,  however, 
present  a  significant  difficulty  in  the  construction  of  the  superharmonic  function  h. 
Can  a  superharmonic  function  h  be  found  that  satisfies 

ft(l)  ■  \\*  -  i^w->mx)) 

for  all  x  E  Rd  -  C,j  =  1,2,3,...?  More  specifically,  can  we  find  a  superharmonic 
function  which  satisfies  the  above  property  without  resorting  to  an  infinite  sum  as  we 
have  to  do  in  this  study?  It  is  unclear  even  in  the  two  ball  case. 

Another  aspect  concerning  the  distance  between  balls  is  the  structure  of  the 
horizon  set  H,  from  a  pure  set  theory  point  of  view.  This  can  probably  be  better 
illustrated  in  R2.  If  three  non-overlapping  open  discs  Bx,  B2  and  B3  all  have  positive 
distance  between  them,  the  set  H,  as  the  result  of  iterated  inversions,  is  the  inter- 
section of  a  sequence  of  nested  closed  sets.  It  is  a  subset  of  a  circle  B  intersecting 
Bi,  B2  and  B3.  The  linear  Lebesgue  measure  of  H,  on  B,  is  zero.  However,  if  each 
pair  of  dBi,  dB2  and  dB3  intersect  at  one  point,  i.e.,  if  these  three  open  discs  are 
"touching"  one  another,  then  the  resultant  H  is  the  entire  circle  passing  through  the 
three  points 

dBx  n  dB2,  dBx  n  dB3  and  dB2  n  dB3. 
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